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Abstract 

We present a new first-principles theory of dynamic arrest in colloidal mixtures based on the 
multi-component self-consistent generalized Langevin equation (SCGLE) theory of colloid dynam- 
ics [Phys. Rev. E 72, 031107 (2005); ibid 76, 039902 (2007)]. We illustrate its application with the 
description of dynamic arrest in two simple model colloidal mixtures, namely, the hard-sphere and 
the repulsive Yukawa binary mixtures. Our results include the observation of the two patterns of 
dynamic arrest, one in which both species become simultaneously arrested, and the other involving 
the sequential arrest of the two species. The latter case gives rise to mixed states in which one 
species is arrested while the other species remains mobile. We also derive the ("fixed point") equa- 
tions for the non-ergodic parameters of the system, which takes the surprisingly simple form of a 
system of coupled equations for the localization length of the particles of each species. The solution 
of this system of equations indicates unambiguously which species is arrested (finite localization 
length) and which species remains ergodic (infinite localization length). As a result, we are able 
to draw the entire ergodic-non-ergodic phase diagram of the binary hard-sphere mixture. 

PACS numbers: 64.70.Pf, 61.20.Gy, 47.57.J- 
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I. INTRODUCTION 



The fundamental understanding of dynamically arrested states is a major challenge of 



contemporary statistical physics and materials science 



2j,|3j. Model experimental colloidal 



suspensions, whose dynamics has been the subject of study on its own right 



have 



played an essential role in the study of dynamic arrest phenomena, providing experimental 



realizations in finely-controlled systems and conditions [7|, [8|, [9|, [10J, [11]. Thus, it is an exper- 
imental fact that by changing a macroscopic control parameter such as the volume fraction, 
or by tuning the effective inter-particle interactions, one may drive the system from the 
region of equilibrium (ergodic) states to the region of dynamically arrested (non-ergodic) 
states. First-principles models and theories that predict and describe these transitions thus 
constitute an essential aspect of the fundamental understanding of these phenomena. The 
mode coupling theory (MCT) of the ideal glass transition j^, 12] is perhaps the most com- 
prehensive theory of this sort, some of whose predictions have found beautiful experimental 



confirmation 



ll| . Many issues, however, still remain to 



)e understood, and this has 



15) . or the development 



14 



prompted the proposal of extended versions of this theory [131 . 
of alternative approaches [if]]. 

In this context, in recent work a new first-principles theory of dynamic arrest has been 



proposed 



17 



18 



19) which is essentially the application of the sel 
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-consistent generalized 



23j to the analysis of 



Langevin equation (SCGLE) theory of colloid dynamics [20l . 
dynamic arrest phenomena in colloidal systems. The SCGLE theory was originally devised to 
describe the tracer and collective diffusion properties of colloidal dispersions in the short- and 
intermediate-times regimes. Its self-consistent character, however, introduces a non-linear 
dynamic feedback, leading to the prediction of dynamic arrest in these systems, similar to 
that exhibited by the mode coupling (MC) theory of the ideal glass transition. The resulting 
theory of dynamic arrest in colloidal dispersions was applied in recent work to describe the 
glass transition in three mono-disperse experimental systems with specific (hard-sphere, 
screened electrostatic, and depletion) inter-particle effective forces Q], with the same 
quantitative predictive power as conventional MCT, but with a smaller degree of difficulty 
in its application 19) . The present paper introduces the multi-component extension of the 



SCGLE theory of colloid dynamics 
dynamic arrest in colloidal mixtures. 



24 



251 ] as a new alternative first-principles theory of 
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Co 
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23, |2g,|22, 



loidal mixtures offer a richer variety of possible dynamic arrest scenarios 
3l| . These include, for example, the possibility that the particles of one species become 
arrested, while the particles of the other species continue to diffuse through the disordered 
matrix of the arrested par ticles. In fact, an important prediction of the multi-component 



extension of MCT 



32 



331 ] was precisely the existence of these mixed states 
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Later studies of dynamic arrest in mixtures [37|, |38|, 



39 



40 



41 



36]. 



42J have been based on MCT 



not because it is the ultimate and definitive theory, but because it was virtually the only 
first-principles theory of dynamic arrest available that had a high level of predictive capa- 
bility. However, relevant issues remain that have not been resolved by its application. For 
example, although the one-component version of MCT predicts the existence of repulsive 
glasses in a hard-sphere system, and its melting upon the addition of a small amount of 
smaller particles (depletion forces), followed by the reentrance to a new kind of (attractive) 
glass 10j, no reentrance seems to be predicted by the two-component version of MCT 43]. 
Issues like this merit the development of alternative approaches to the description of dynamic 
arrest phenomena in mixtures, particularly if they are independent from conventional MCT. 
Although the SCGLE theory presented here has several features in common with MCT, 
such as its use of exact memory function expressions for the concentration-concentration 
time-correlation functions, in the SCGLE theory the derivations, approximations and ar- 
gumentation are conceived and based on a totally different fundamental and conceptual 
framework 
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23 



24 



25]. 



Besides presenting the SCGLE theory of dynamic arrest in colloid mixtures, in the present 
paper we illustrate the scenarios predicted by this theory in the simplest conditions, namely, 
in its application to model binary mixtures. For two such systems, a repulsive Yukawa 
mixture and a hard-sphere mixture, we numerically solve the full self consistent system 
of dynamic equations for the so-called propagators (or Green's functions) of the self and 
collective intermediate scattering functions. We illustrate in this manner the two possible 
dynamic arrest patterns possible in a binary colloidal mixture, namely, the simultaneous 
and the sequential arrest of the two species. In the latter case, we discuss the hybrid or 
partially-arrested states in which one species is dynamically arrested while the other species 
remains ergo die. 

Determining the location of the ergodic-non-ergodic transitions and calculating the non- 
ergodic parameters may be done by solving the self-consistent system of equation or, more 
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simply, by finding the long-time asymptotic stationary ("fixed-point") solutions of the self- 
consistent system of equations. Here we also address this issue, and report the derivation of a 
strikingly simple and general result which takes the form of a system of coupled equations for 
the localization length of the particles of each species. Solving this system of equations leads 
to the unambiguous identification of which species is arrested (finite localization length) and 
which species remains ergodic (infinite localization length). As a result, we are able to draw 
the entire ergodic-non-ergodic phase diagram of the binary hard-sphere mixture, which is 
one of the main specific contributions of this work. 

The present paper starts with a brief summary of the SCGLE theory of the dynamics 
of colloidal mixtures. This theory was developed in Ref. 25| and tested there through its 
comparison with the short- and intermediate-time results of Brownian dynamics computer 
simulations for a model colloidal mixture interacting through screened Coulomb pair inter- 
actions. That original version of the SCGLE theory required as an input not only the partial 
static structure factors, but also other equilibrium static structural information associated 
with the exact short-time conditions that this initial version of the theory aimed at incorpo- 
rating. As it has been shown more recently [l9], however, a simplified version of the SCGLE 
theory, in which this exact short-time information is neglected, happens to be equally accu- 
rate, particularly in the long-time regime and with respect to the predicted dynamic arrest 
scenarios. Thus, we shall refer to this simplified version simply as the SCGLE theory, which 
we review in the following section. 

Although our main interest here is in the description of dynamic arrest phenomena, in 
section III we start our discussion by reviewing the predictive capability of this simplified 
theory in the short- and intermediate-time regimes of the model binary Yukawa mixture 
studied in Ref. in one of the fully ergodic states considered in that work. We then 
proceed to illustrate the dynamic arrest scenarios predicted by this theory, by taking as 
the initial reference the same repulsive Yukawa mixture, and increase the intensity of the 
repulsion parameter (a process that amounts to effectively lowering the temperature), until 
reaching a transition to a dynamically arrested state. This turns out to be a transition from 
a fully ergodic to a fully arrested state in which both species are simultaneously arrested. 
In the second of our illustrative examples, in section IV we consider a hard-sphere model 
colloidal mixture, whose state parameters are the size disparity and the concentration of 
the two species. Here we present an example of a transition to a partially arrested state, 
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in which the large particles become arrested and the smaller particles continue to diffuse in 
the disordered matrix formed by the arrested species. In section V we derive the equations 
for the non-ergodic parameters of the system, and illustrate its use in the context of the 
hard-sphere binary mixture by determining its entire dynamic-arrest phase diagram. The 
main conclusions of this paper are finally summarized in section VI. 



II. THE SCGLE THEORY OF THE DYNAMICS OF COLLOIDAL MIXTURES 



In this section we summarize the main elements of the self- consistent theory of the 
dynamics of colloidal mixtures presented in Refs. 



24 



For a colloidal mixture with v 



species, the dynamic properties can be described in terms of the relaxation of the fluctuations 
Sn a (r, t) of the local concentration n a (r, t) of colloidal particles of species a (= 1,2, ...v) 
around its bulk equilibrium value n a . The average decay of Sn a (r, t) is described by the time- 
dependent correlation function F a p(k,t) = (5n a (k, t)5rip(— k, 0)) of the Fourier transform 
5n a (k, t) = (l/y/N a ) YliLi ex P ' r | a ' > (^)] °f the fluctuations 8n a (r, t), with r[ a \t) being 
the position of particle i of species a at time t. F a/3 (k,t) is referred to as the partial 
intermediate scattering function, measured by experimental techniques such as dynamic 
light scattering 4j, |6[]. One can also define the self component of F a p(k,t), referred to 
as the self-intermediate scattering function, as F$j(k,t) = 5 aP (exp [ik • AR (a) (*)]), where 
AR^'(t) is the displacement of any of the N a particles of species a over a time t, and 5 a p 
is Kronecker's delta function. 

For notational economy, let us denote the matrix with elements F a p(k,t) (at, ft = 
1,2,.. .,v) by F(k,t) or simply by F(t), and similarly for F^ s \t). These matrices have 
the initial values F(t = 0) = S and F^ s \t = 0) = I, where I is the identity matrix, 
I a p = d^, and the elements of the matrix S are the partial static structure factors, 
S Q p(k) = (8n a (\t, 0)5rip(— k, 0)). The matrix of intermediate scattering functions may also 
be normalized by its initial value, to define the collective propagator matrix = 
which is such that ^(0) = /. In a similar manner, we define the self propagator matrix as 
*W(t) = FW(t), normalized to tfW(O) = /. 

The multi-component version of the SCGLE theory starts with the exact time-evolution 
equation that governs the relaxation of the Fourier transform (FT) 5n a (k,t) of the fluc- 
tuations of the local concentration of species a. This equation also governs the relaxation 
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of the partial intermediate scattering functions F a p(k,t), and can be written as an exact 
expression for the Laplace transform (LT) F(z) of the matrix F(t), namely [241 ]. 

F(z) = {z + (I + C{z)Y 1 k 2 DS- l Y l S, (2.1) 

where the elements of the matrix D are given by D a p = SapD^, with being the diffusion 
coefficient of species a in the absence of interactions. This is related with the solvent 
friction coefficient on an isolated particle of species a, through the Einstein relation, 

the matrix C(z) are the LT of the so-called 



44|. The corresponding result for the "self" 



= ksT/C^. The elements C a /3(k,z) o 
irreducible memory functions C Q/ g(/c,t) [36 
component, F^(z), of the matrix F(z) reads 

= {z + (I + C^ s \z))- l k 2 Dy l , (2.2) 

where the matrix C^ B \z) is the corresponding irreducible memory function. 

The second ingredient of the SCGLE theory is the intuitive notion that collective and self 
dynamics may be connected in a simple manner. Vineyard's approximation, in which \I/(t) 
is approximated by \E f ^- ) (t), is a rather primitive implementation of this idea. The SCGLE 
theory incorporates the same notion, but implemented at the level of the irreducible memory 
functions C(t) and C^(t), namely, 

C(t) = C (s) (t), (2.3) 
which we shall simply refer as the Vineyard-like approximation. Let us mention that in the 



original proposal of the mult i- component version of the SCGLE theory 25[ the difference 
[C(t) — C^(t)] was approximated not by zero, as in this equation, but by its exact ("single 
exponential" ) short-time limit [C SEXP (t) —C^' SEXP (t)}, given the fact that explicit expres- 
sions exist 0, Q] for the short-time limits C SEXP (t) and C^' SEXP (t) of the irreducible 
memory functions C(t) and C^ s \t) in terms of equilibrium static structural properties of 
the system, which are considered known. More recently, however, it was found [19( that an 
equally accurate but much simpler approximation results if the difference [C(t) — C^ s '(t)] is 
approximated by zero, as in Eq. (I2.3p . thus eliminating the need to determine the short-time 
memory functions C SEXP {t) and C i - s) ' SEXP {t). 



The third ingredient of the SCGLE theory was suggested |25j as the interpolation of 
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C( s '(t) between its two exact limits at small and large wave- vectors, by means of a phe- 
nomenological but universal interpolation device. We know that the C^> SEXP (t) is also 
the exact large-fc limit of C^ s \t), whereas the small- A; limit is the so-called time-dependent 
friction function matrix AC*(t). Thus, we write C {s \k,t) = C (s) > SEXP (k,t) + [A(*(t) - 
C^ s '' SEXP (k,t)]\(k), where we have written explicitly the ^-dependence of the various ma- 
trices. The matrix A is a matrix of phenomenological interpolating functions, whose elements 



are defined as 



22 



25] 



A a/ ^ = W[l + (VCiJ 2 ], (2-4) 

with h^ n being the position of the first minimum following immediately the main peak of 
the partial static structure factor S aa (k). At times long enough that the short-time memory 
function C^ ,SEXP (k,t) has already decayed, the interpolation formula above for C s (k,t) 
simplifies still further, reading 



C^(k,t) = [AC(t)]X(k). 



(2.5) 



As suggested in Ref. [19(, it is this approximation that enters in the present formulation of 
the multi-component SCGLE theory of dynamic arrest. 

The matrix A(*(t) in eq. ( 12.5ft is a diagonal matrix whose oeth diagonal element, A£*(£), 
is the time-dependent friction function of particles of species a. The final ingredient of the 
SCGLE theory is an approximate but general result for this property, which, as explained 



below, was derived as an application of the GLE formalism 
phenomena. Snch an expression reads f 



46 



43, 



481 ] to tracer diffusion 



Di 



ci 3 fcA; 2 [F (s) (t)] Qa [cv / ^F(t)5- 1 v / ^^ c 



(2.6) 



3(2tt) 3 

with the elements of the ^-dependent matrices h and c being the Fourier transforms h a p(k) 
and c a p{k) of the Ornstein-Zernike total and direct correlation functions, respectively. Thus, 
h and c are related to S by S = I + y/nh^/n — [I — ^Jncy/n\~ 1 , with the matrix y/n defined 
as [y/n\ a /3 = S a/3 y/n^. 



Let us now comment on the origin of Eq. ( 12.6H . In Ref. 46j the effective Langevin 
equation of a tracer colloidal particle interacting with the other particles of a mono-disperse 
suspension was derived, using the concept of contraction of the description {47] (a summary 
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of such derivation is contained in Appendix B of Ref. 18] ). Besides the solvent friction 
force, — C°V(i), the direct interactions of the tracer particle with the other particles give rise 
to an additional friction term of the form — J * dt'A((t — t')V(t'), where V(t) is the tracer 
particle's velocity at time t. In the process, an exact result for the time-dependent friction 
function A£*(i) = A((t)/(° is generated which, when complemented with two well defined 
approximations, leads to the mono-disperse version of Eq. (12.61) . The multi-component 



extension was derived in Ref. 48J] from the exact result in eq. (124) of that reference, upon 
the introduction of the "homogeneous fluid" and the "decoupling" approximations (eqs. 
(135)- (138) and eq. (144) of the same reference). 

Finally, let us mention that eqs. (12.1 p and (12. 2p are exact results, and that eq. (12. 6p 
derives from the exact result just referred to. Hence, it should not be a surprise that the 
same results are employed by other theories; in fact, these equations also appear in the 
formulation of MCT. The difference lies, of course, in the way we relate them and use them; 
in this sense, the distinctive elements of our theory are then the Vineyard-like approximation 
in eq. ( 12.31) and the factorization approximation in eq. ( 12.51) . 

III. ILLUSTRATIVE APPLICATION: THE MODEL YUKAWA MIXTURE 

The set of coupled equations ( 12. II) - ( 121)1) constitutes the self-consistent generalized 
Langevin equation theory of the dynamics of colloidal mixtures, and has to be solved nu- 
merically. For this, equations (12. ip and (12. 2p are first Laplace-inverted, and written as a set 
of coupled integro-differential equations involving functions of k and t. These functions are 
then discretized in a mesh of points large enough to ensure independence of the solution 
with respect to the size of the mesh. The discretized system of equations is then solved 
by a straightforward direct iteration method or, for the long times req uired in the study of 



25] to 



dynamic arrest phenomena, by the methods described in Refs. [49|, [50 ] . 

To illustrate its quantitative accuracy, the SCGLE theory was applied in Ref. 
the calculation of the dynamic properties of a particular model system, namely, a binary 
mixture of Brownian particles interacting through a hard-core pair potential of diameter a, 
assumed the same for both species, plus a repulsive Yukawa tail of the form (in units of the 
thermal energy ksT = /3 _1 ) 
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f3u aP {r) = A a Aj3 -j- . (3.1) 

The dimensionless parameters that define the thermodynamic state of this system are the 
total volume fraction = |n<7 3 (with n being the total number concentration, n = n\ + ^2), 
the relative concentrations x a = (a = 1,2), and the potential parameters Ai, A 2 , and 
z. As in Ref. 25[], here we also assume that the free-diffusion coefficient of both species 



are identical, i.e., D\ = D® = D°. Explicit values of the parameters a and D° are not 
needed, since the dimensionless dynamic properties, such as F a/3 (k,t), only depend on the 
dimensionless parameters specified above, when expressed in terms of the scaled variables ha, 



and t/to, where t = a 2 /D°. In Ref. 25j Brownian dynamics simulations were reported for 
the static and dynamic properties of this system with fixed screening parameter z = 0.15 and 
various coupling parameters A\ and A2 and volume fractions, in the short- and intermediate- 
time regimes. 

The partial static structure factors S a p(k) obtained from those simulations are also em- 
ployed here as the static input needed in the implementation of the SCGLE theory. As a 
result, we can calculate all the dynamic properties entering in the self-consistent equations 
above, and other properties that derive from them, such as the mean squared displacements 
W a (t) =< [AR^(t)] 2 > /6 or the time- dependent diffusion coefficients D a (t) = W a {t)/t. 
In Fig. [1] a comparison of these results with the Brownian dynamics calculations is pre- 
sented, involving the case of a binary mixture with z = 0.15, A\ = 10, and A 2 = lOy^- 
The volume fraction of the more interacting species is kept fixed at 02 = 2.2 x 10 -4 and 
0i takes the values 0i = 0.725 x 10 -4 (left column), 0i = 2.2 x 10~ 4 (middle column), and 



6.6 x 10 (right column), corresponding to molar fractions x% = 0.25, X\ = 0.5 and 



x\ = 0.75, respectively. This figure corresponds to the same conditions as Fig. 3 of Ref. [25], 
and illustrates the fact that the present simplified SCGLE theory also provides an excellent 
description of the relaxation of concentration fluctuations in colloidal suspensions in the 
short- and intermediate-time regimes illustrated in the figure. We notice that the largest 
departures of the theoretical results from the simulation data appear only at the short times 
illustrated by the time t = to, but the agreement improves at longer times, as illustrated 
by the time t = 10to- The main message of this figure is that the present self-consistent 
theory of colloid dynamics provides an adequate description of the dynamic properties of a 
simple idealized mixture also in the short- and intermediate-time regimes, even though it 



does not satisfy the exact short-time conditions built in the original proposal of this theory 



25j | . With this confidence, let us now explore the long-time regime under conditions typical 



of the vicinity of a transition to dynamically arrested states. 

For this, let us take the conditions of the middle column of Fig. 1 as a reference, and let 
us keep the parameters 0i, <p 2 and z fixed (0i = <fi 2 = 2.2 x 10 -4 and z = 0.15) while we 
increase the coupling parameters A\ and A 2 in an identical proportion, to mimic a process 
of lowering the temperature. Thus, let us write A a = A° a x (a = 1,2), with A\ = 10 and 
A° 2 = 10^5, and increase x from its unit value in Fig. 1 until a transition of dynamic 
arrest is encountered. In this case, we do not run computer simulations to provide the input 
partial static structure factors needed in the self-consistent scheme. Instead, we resort to 
an approximate liquid state theory, namely, the solution of the Ornstein-Zernike integral 
equation within the hyper-netted-chain (HNC) approximation 511 ] . In doing this we only 
loose quantitative precision, but no qualitative accuracy. In fact, the HNC approximation 
will yield the same scenario, but at effective values of the coupling parameters Ai and A 2 
larger than the actual ones. For example, the fitting curves for the partial static structure 
factors in the middle column of Fig. 1 were obtained using the HNC approximation with 
effective coupling parameters given by A a = A Q a x with x — 1-225. 

In Fig. [2] we present the results of the SCGLE theory for the diagonal elements ty aa (k, t) of 
the collective propagator matrix \I/(t) as a function of time for fixed wave- vector k = k max , 
where k max is the position of the main peak of the corresponding static structure factor 
S aa (k). This figure illustrates the fact that as we lower the temperature starting from the 
reference state (x = 1.225) the system reaches eventually a transition to an arrested state. 
This occurs for Ai = 42.57, and A 2 = 95.19, i.e., for x — 4.257, and this is a transition from 
a fully ergodic to a fully arrested state, i.e., both species become simultaneously arrested. 
This is evidenced by the fact that below this threshold condition, exemplified in the figure 
by the case x = 4.255, both propagators decay to zero, whereas above this condition, 
exemplified by the case x = 4.330, they decay to a finite non-ergodic value. Thus, the non- 
ergodic parameters jump discontinuously at the transition, from zero in the ergodic side, to 
finite values in the glass side. Exactly the same behavior is exhibited by the self-diffusion 
propagators £) and W^i^i 0> no ^ included in the present illustration. Let us stress 

that this process of increasing A\ and A 2 keeping z fixed is a mathematical, rather than 
physical, exercise; the practicable experimental manner to increase A\ and A 2 would be to 
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increase the particle's electric charge (or to lower the dielectric constant of the solvent), 
which in reality would also increase the screening parameter z. Here, however, we use this 
exercise only to illustrate the pattern of dynamic arrest transition in colloidal mixtures in 
which the particles of both species become localized simultaneously. 

Similar calculations can be carried out for other regions of the parameter space 
(<f)i,(f)2,Ai,A2,z), but this scanning and its results in terms of other possible regions that 
exhibit interesting dynamic behavior constitutes an independent issue, beyond our current 
objective of introducing our SCGLE theory of dynamic arrest in mixtures and illustrating its 
use. We can mention, however, that in other regions of this state space the theory predicts 
also the other pattern of dynamic arrest, which we now illustrate in the context of a simpler 
system, namely, the hard-sphere colloidal mixture. 



IV. THE HARD-SPHERE MIXTURE 



Let us now consider the binary hard spheres model dispersion, i.e., the uncharged version 
of the previous system in which we now allow for hard-sphere size disparity. Thus, we 
consider two species of particles of hard core diameters o\ and <r 2 , present at concentrations 
rii and n 2 . We use as dimensionless control parameters the size asymmetry 5 = (J\j(Ji < 1 
and either the two volume fractions 4>i an d 02 (with Q = iin a a^/6), or the total packing 
fraction = 0i+02 together with the molar fraction of the smaller species, x\ = 711/(711+712)- 
This system serves to illustrate some features associated precisely with the size disparity. 
As discussed in more detail below, one expects that for mild asymmetries the only possible 
mode of dynamic arrest for all volume fractions and all molar fractions is that in which 
both species are arrested simultaneously. It is only in the regime of disparate sizes that one 
expects a richer scenario, including the second mode of dynamic arrest, in which the large 
particles are first arrested while the other species continues to diffuse, followed by a second 
transition at a higher total volume fraction in which the smaller particles are finally localized 
inside the pores of the matrix formed by the previously arrested particles. Thus, there must 
be a threshold asymmetry 5 C beyond which this might occur. One also expects, however, 
that even in the regime of considerable size disparities but for conditions corresponding to a 
few large particles in a see of smaller ones, the only mode of dynamic arrest continues to be 
the mode of simultaneous dynamic arrest of the two species, since in this case the arrest of 
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the small particles implies the arrest of the few large ones. We call this the "chancaquilla' 



limit 



52( | . in reference to one of the many names of a rather universal sweet prepared by 



the vitrification of liquid molasses in which nuts or raisins are previously dispersed. Clearly, 
the dynamic arrest of the supporting molasses implies the arrest of the nuts dispersed in it. 
Hence, for large size asymmetries we expect to find both modes of dynamic arrest located 
in the neighborhood of the opposite composition limits, x± = and x\ = 1. 

In fact, this is exactly the scenario predicted by our theory. We illustrate this with the 
SCGLE results obtained using the Percus-Yevick [53J structure factors as the static input, for 
a size asymmetry 5 = 0.3. Let us calculate the collective and self propagator matrices ^f(t) 
and #W(t) at increasing total volume fractions starting in the ergodic regime, crossing the 
dynamic arrest transition line, and ending inside the regime of arrested states. If this is done 
in the neighborhood of the X\ — 1 axis, say for x\ = 0.8, the scenario would be essentially 
that described by the chancaquilla mechanism involving the simultaneous dynamic arrest of 
both species, with similar results to those discussed in Fig. [2] for the Yukawa mixture. The 
other mode of dynamic arrest is observed in the opposite regime, where the small particles 
are a minority, exemplified by the molar fraction x\ = 0.2. This case is illustrated in Fig. 
|3l which presents the evolution of the relaxation of the diagonal elements of the collective 
and self propagator matrices *ff(t) and ty( s \t) as we increase the total volume fraction 
at fixed molar fraction X\ = 0.2, starting in the ergodic regime. We find that, rather than 
passing directly from the fully ergodic region to the region of fully arrested states, as in Fig. 
El one first passes through an intermediate region of hybrid states, in which only the big 
particles are arrested but the small particles continue to diffuse. Thus, at a volume fraction 
(f)y £s 0.545, the system passes from the fully ergodic region to the region (j)g < <p < (ffi 

(2) 

of hybrid states, with <\> g being the location of the second transition consisting of the 
localization of the small particles in the pores of the random matrix of big particles. In the 
present example this occurs at 4%' « 0.555. 

Let us mention that in all these cases one observes that the diagonal collective diffu- 
sion propagator \P aa (fc, £) behave in qualitatively identical manner to the corresponding self 
diffusion propagator Waa(k,t) in the sense that either both relax to zero or both exhibit 
dynamic arrest (i.e., decay to a finite value). One does not encounter a situation in which, 
for example, ^ ii ( A; , t) decays to a finite value while ty^(k,t) decays to zero. This indicates 
that at the level of the diagonal propagators, self and collective diffusion provide consistent 
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descriptions of the pattern of relaxation. Thus, we may identify the fully ergodic states by 
the fact that all the elements of the (self and collective) propagators decay to zero, whereas 
the fully non-ergodic states can be identified by the decay of all the propagators to finite 
asymptotic values (identified with the non-ergodic parameters ip and ip s ). The results for 
= 0.53 in Fig. [3] are representative of fully ergodic states, whereas those for = 0.57 il- 
lustrate fully arrested states. Hybrid states are characterized by a mixed behavior, in which 
the diagonal propagators ^n(k, t) and ty$(k,t) associated to the mobile species, decay to 
zero, while the diagonal propagators ^22(^5 1) and ^f^(k,t), associated to the large par- 
ticles, exhibit arrested behavior. These conditions are illustrated in Fig. [3] by the results 
corresponding to = 0.55. 

The relaxation of the crossed collective propagators ^12 (/c, t) and ^2i{k,t), not shown 
in Fig. [31 do deserve some discussion, particularly in connection with mixed states. First, 
one must have in mind that, although the matrices F(k,t) and S(k) are symmetric, their 
product F(k, t)S~ 1 (k) = *$?(k, t) in general is not; thus, in general, ^12 (A;, t) 7^ ^2i(k,t). 
Second, since the initial value ^(k, 0) of the propagator matrix is the unit matrix, the initial 
value of the off-diagonal collective propagators ^12 (A;, i) and ^2i(k,t) must vanish. Third, 
the off-diagonal collective propagators are not necessarily positive. This is illustrated in 
Fig. HI which exhibits the time-relaxation of the four collective propagators ty a p(k,t) for 
the same mixed state included in Fig. [3l but now evaluated at a common wave-vector k. 
The wave-vector employed in this figure corresponds to the position of the maximum of 
S22(k); for reference, the inset of this figure contains the partial static structure factors of 
this system. At this wave-vector we observe that both, *&i2(k, t) and \&2i(&, t), are negative, 
and that *$?i2(k,t) relaxes to a finite non-ergodic asymptotic value, while ^2i{k,t) is always 
much smaller and relaxes to zero. 

Let us mention that one could also display the same information not in terms of 
propagators but in terms of the so-called collective "correlators" $ a p(k,t), defined as 
&ap(k,t) = F a p(k,t)/S a p(k). Thus defined, the diagonal collective correlator & aa (k,t) 
is a linear combination of the corresponding diagonal collective propagator ty aa (k,t) and 
of a crossed propagator; more precisely, $n(/c,t) = ^u(k, t) + [S2i(k)/ Sn(k)]^/ 12 (k, t) and 
$22(^,0 = ^22(^5 1) + [Si2(k) / S22(k)]^2i(k, t). In the latter case, since \& 2 i(fc, £) is small 
and relaxes to zero, we find that the collective correlator $22^, t) is virtually identical to 
the collective propagator ^22^, t), as shown in Fig. HI In contrast, although the propagator 
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(k, t) does decay to zero, the diagonal collective correlator $n(fc,t) exhibits dynamic 
arrest due to its linear dependence on ^>i 2 (A;,t) which, as illustrated in the figure, relaxes 
to a non-zero value. This observation may be relevant when describing the hybrid states 
above. For example, the description in terms of correlators has been employed in recent 
work 30, |3g, |42|], where the situation illustrated in Fig. H] arises, leading to the notion 



that collective diffusion in mixtures is intrinsically different from self-diffusion. Our results 
indicate, however, that this apparent difference is just the result of the normalization con- 
vention employed in the definition of the correlators, and that if the same information were 
displayed in terms of the elements of the propagator matrices, no mixing of relaxation modes 
would occur that suggest this notion. We should stress that the self correlators, defined 
as $^l(k,t) = F^g(k,t), are identical to the self propagators ^l(k, t), and hence correctly 
describe the fact that only the small particles diffuse and the large particles are arrested. 

Finally, let us mention that distinguishing the three dynamic states above (fully ergodic, 
fully non-ergodic, and hybrid states) with only the information of the type in Figs. [3]and|4] 
is not a simple matter. The reason is that the non-ergodic parameters of the propagators 
associated with the localization of the small spheres do not jump discontinuously from zero 

(2) 

to a finite appreciable value at <pg . Instead, they increase continuously from a value of 
zero right at this transition, and hence, attain finite but very small values in the close 
neighborhood of the transition. In the following section, however, we shall derive a simpler 
procedure to precisely locate these transitions and to evaluate the non-ergodic parameters 
in the (partially or totally) arrested states. 



V. FIXED-POINT EQUATIONS 

With the information just presented, one may choose to analyze the detailed time- 
dependence of the propagators, and discuss the various regimes of the relaxation processes 
described by these functions at specific points in the state space. Alternatively, one can 
discuss more global properties, such as the general topology of the ergodic-nonergodic phase 
diagram of the system, just like one does in equilibrium statistical thermodynamics. In fact, 
the ergodic-nonergodic phase diagram should be an important complement to conventional 
equilibrium phase diagrams in terms of practical applications. If one is interested specifi- 
cally in such global properties, in principle one does not need to solve the set of coupled 
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equations (I2.1l) - ()2.6p for the full time- and wave- vector-dependence of all the dynamic prop- 
erties involved. Instead, one may solve only the so-called fixed-point equations, i.e., the 
equations for the long-time asymptotic stationary solutions of the full SCGLE theory. The 
solutions thus determined are precisely the non-ergodic parameters, whose zero or non-zero 
value reveal if the state of the system is ergodic or dynamically arrested. 

In order to derive the fixed point equations, let us define the non-ergodic parameter 
of the dynamic properties ^(t), ^ s (t), C(t), C s (t), and A(*(t), as the matrices if)(k) = 
linw^M), r(k) = linw^tM), c(k) = linw C(k, t), c s {k) = lim^ C s (k, t), 
and A^** 00 ** ee lim^oo A£* (t) , respectively. All the elements of these matrices are zero in a 
fully ergodic state, and the non-zero value of some or all of them indicate partial or total 
non-ergodicity. Let us now substitute the dynamic properties involved in the self-consistent 
system of equations in Eqs. (I2.1l) - ()2.6p by the sum of their asymptotic long-time value 
above, plus the rest (which, by definition, always relaxes to zero). In the resulting system 
of equations, let us then take the asymptotic long-time limit, thus generating the following 
self-consistent system of equations for the non-ergodic parameters 

■ip(k) = [c(Jfe) + tfDS- 1 ^)}- 1 ^), (5.1) 

*p s (k) = [c s (k) + k 2 D}- l c s (k), (5.2) 

c{k) = c s {k), (5.3) 

c s (k) = \{k)AC (QO \ (5.4) 

and 

AC (oo) = J d 3 kk 2 [^% a [c^V^h} aa , (5.5) 

where the matrices c and h inside the integral in the last equation refer to the direct and total 
Ornstein-Zernike correlation functions, and together with the matrix ^/n, were defined below 
Eq. (12.6j) . Using Eqs. (I5.3P and (15 Ah in Eqs. ( 15.11) and ( 15.21) . we can express the non-ergodic 
parameters ip(k) and ips(k) in terms of A^*(°°\ Substituting the resulting expressions in 
Eq. (15.51) . we finally get the following closed equation for AC*(°°) 
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AC 



*(oo) 
a 



D° 



3(2tt) ; 



d s kk 2 



(A(A;) + A; 2 (AC (oo) )" 1 ^) 1 X(k) 



c\/n 



(A(Jfe) + fc 2 (AC* (oo) )~ 1 J DS'- 1 (jfe)) 1 A(fc)l y/Eh 



(5.6) 



for a = 1, 2, v, with v being the number of species. Clearly, this equation always admits 
the trivial solutions A(a^ = 0, which corresponds to the fully ergodic fluid state. The 
existence of other non-zero real solution(s) is associated to the existence of arrested states. 

It is convenient to re-write Eq. (15. 6p in terms of the diagonal matrix 7 = (A^ 00 )) -1 D , 
whose diagonal elements 7 a will then satisfy the following set of equations 



1 

la 



J d^^I + k^X^ik)] ^^[cV^il + k^X-^S-'ik)] l V^h 



(5.7) 

For the same reasons given in the mono- component case |17|, [18] , the new order parameters 
7 Q may be identified with the long-time asymptotic value of the mean squared displacement 

(msd) of particles of species a. Thus, 7^ is infinite if that species diffuses, and is finite if it 

1/2 

is arrested; in fact, in the latter case 7 Q is the localization length of the particles of that 
arrested species. In terms of 7, Eqs. ( 15. II) and ( 15.20 for the non-ergodic parameters ip(k) 
and ip s {k) read 



and 



^j(k) = [I + k 2 1 X- 1 {k)S-\k)}~ 1 



(5.8) 



ip s (k) = [I + k 2 1 \-\k)] 1 



(5.9) 



Eqs. ( 15.7p - 115.9p are the fixed point equations of the SCGLE theory. In the particular case 
of a mono-disperse system, one recovers the result derived and employed in Refs. Q] • 
For a given system, i.e., for given pair potentials, one is supposed to determine first the 
static structure factors S a p(k), which is the only external input in Eq. (J57TJ) (recall that the 
matrix X(k) is determined by these structural properties, according to Eq. (I2.4p ). Eqs. (15. 7p 
can be solved numerically, and the solutions for 7 are then employed in Eqs. ( 15.81) and ( 15. 9p 
to determine the non-ergodic parameters. 
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Let us apply Eqs. (15. 7\i to calculate the order parameters 71 and 72 for binary hard-sphere 
mixtures. The corresponding partial static structure factors S a p{k) needed in that equation 
will again be provided by the Percus-Yevick (PY) approximation. From the resulting values 
of 71 and 72, we can classify the states of this system as ergodic, fully non-ergodic, or partially 
non-ergodic, and thus scan the state space to determine the regions of these dynamic phases 
and the boundaries between them. For this, let us adopt the asymmetry parameter S, the 
total volume fraction 0, and the molar fraction x\ of the smaller spheres as the dimensionless 
state variables spanning the state space (5, 0, X\). 

Before presenting the results of this exercise, that we collect in Fig. El let us say that 
from intuitive arguments one may advance some of the features of this phase diagram. For 
example, as already indicated in section IV one expects that for mild asymmetries the only 
possible mode of dynamic arrest for all volume fractions and all molar fractions is that in 
which both species are arrested simultaneously. The reason for this expectation is that this 
regime must contain the extreme case corresponding to the degenerate limit 6 = 1. From 
the point of view of their interactions, in this limit the particles of both species are actually 
indistinguishable (<Ti = 02), thus being in fact a mono-disperse system, whose particles have 
been classified only artificially as belonging to either of these two species. We expect, then, 
that for any molar fraction < X\ < 1, the system will be arrested at the total volume 
fraction (j)^ at which a mono-disperse system would be arrested. This means that in the 
plane (5 = 1, 0, xi), the transition line must correspond to the condition 0i + 02 = 5 = 0g"^ ■ 
According to the mono- component version of the SCGLE theory, and within the same level 
of approximation adopted here for the static structure factor (i.e., the PY approximation) 



we know that = 0.537 



181 ]. Thus, in the degenerate case the straight line 



, (m) 



(horizontal line in Fig. El a) separates the ergodic region < <f)^ from the region > (j)^ 
of totally non-ergodic states, in which both species are dynamically arrested. 

Thus, the first question is now if the SCGLE theory has built-in this degenerate limit, 
i.e., if under these degenerate conditions, Eqs. (I5.7p -( l579l reduce to the mono-disperse result 



discussed in references 
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181 ]. Indeed, if one constrains the solution for 71 and 72 to 



exhibit this symmetry, i.e., that 71 = 72, one can show that the two equations involved in 
Eq. ( 15. 7p are indeed identical to each other and to the corresponding mono-disperse equation, 
independently of the composition represented by the molar fraction x\. This means that 
both "species" either remain ergodic or are arrested together. 
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The next question is then what happens when this degeneration is broken, so that u\ < 02, 
but remaining in the region of mild asymmetry 5 ^ 1, in which we expect small deviations 
from the degenerate limit. Unfortunately, we do not have good intuitive arguments to 
anticipate even the sign of these deviations. In Fig. 0a, however, we report the predictions 
of our theory for this regime, in terms of the glass transition lines corresponding to the 
values of the asymmetry parameter 5 = 1.0, 0.8, 0.6 and 0.4, for which our theory predicts 
that the glass transition line moves to higher total volume fractions <j) g as the size disparity 
increases (i.e., S decreases), except at the two ends X\ = and x\ = 1 of the transition line, 
where the value of <p g must be fixed by its mono-disperse value cj)^ 1 for all size disparities. 

The general trends illustrated by Fig. [5Ja were actually first observed and described by 



Williams and van Megen 



29j in their experimental process of melting an originally mono- 



disperse glass by the replacement of a fraction of its particles by particles of a smaller size 
keeping the same total volume fraction. The dots and the arrow in Fig. [5ja correspond to 
their experimental process, that may conceptually be described as driving a system from 
the initially mono-disperse glass state (5, <p/<f)g m \ x\) = (0.6, 1.0175, 0.0) through the 
arrested states (0.6, 1.0175, 0.196) and (0.6, 1.0175, 0.340), to end in the ergodic state 
(0.6, 1.0175, 0.537) which lies in the ergodic region according to our predictions for the 
same size asymmetry 5 = 0.6 of the experimental system (the representation of this process 
in the scale of Fig. [5j a involves the re-scaling of the experimental volume fractions <p exp = 
0.58 with the experimental value of <% which, within the experimental errors, we took as 
4 m) = 0.57). 

The regime illustrated in Fig. [5] a was also studied by Gotze and Voigtmann with MCT 
371 . l38l |. with predictions qualitatively similar to ours for intermediate size-disparities (8 < 



0.65), but with conflicting predictions for milder size disparities (5 ~ 0.8), where MCT 
predicts "S-shaped" transition lines that may go below the mono-disperse glass transition 
volume fraction 37]. The SCGLE theory, predicts, instead, that the transition lines 
for the mixture, as illustrated in Fig. Oa always lie above the horizontal transition line 
corresponding to the degenerate mono-disperse case. These conflicting predictions of MCT 
and the SCGLE theory could be resolved by experiments and/or computer simulations. The 
actual magnitude of the difference between the predictions of the two theories is, however, so 
small that a conclusive discrimination between them will require experimental or simulation 
results of higher resolution than provided by the currently available data 0, Q] . 
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Fig. [5jb illustrates the regime of large size disparities. The SCGLE theory predicts that 
below a threshold asymmetry 5 C ~ 0.4, a region of mixed states appear, which is illustrated 
by the shaded area in the results for 5 = 0.3. This region is bounded on the left by the 
axis X\ = 0, below by the transition of dynamic arrest of the large particles (solid line), and 
above by the localization transition of the small particles (dashed line). The region of mixed 
states ends on the right in the bifurcation point where these two transition lines meet, and 
where the transition of simultaneous dynamic arrest of both species starts, ending in the 
axis X\ — 1 (solid line). The area of the shaded region increases as 8 decreases, as illustrated 
by the results for 5 = 0.2, in which the region of mixed states was not shaded. Of course, as 
S approaches the threshold value 5 C ~ 0.4, the bifurcation point moves to the axis x\ = 0, 
and the region of mixed states disappears. Although we cannot in practice explore the limit 
5 « 1, we do not expect additional qualitatively different trends beside those exhibited by 
the illustrative cases in Fig. [5j Thus, to our knowledge, this is the first time that the entire 
ergodic-non-ergodic phase diagram of a binary colloidal mixture has been outlined. 



VI. CONCLUDING REMARKS 



In this paper we have introduced the self-consistent generalized Langevin equation theory 
of dynamic arrest in colloidal mixtures, and we have applied it to the description of dynamic 
arrest phenomena in two simple model systems, namely, a repulsive Yukawa and a hard- 
sphere mixtures, for which we illustrated the scenarios predicted by this theory through the 
full numerical solution of the coupled system of dynamic equations for the propagators \I/(t) 
and \l/^(i) of the self and collective intermediate scattering functions. In this manner we 
identified and illustrated the two possible patterns of dynamic arrest in a binary colloidal 
mixture, namely, the simultaneous and the sequential arrest of the two species. The for- 
mer corresponds to the transition from the region of fully ergodic states, characterized by 
infinite values of the order parameters 71 and 72, to the region of fully non-ergodic states, 
characterized by finite values of these two parameters. This transition is characterized by 
a discontinuous jump of both order parameters. This is the only dynamic arrest transition 
possible for moderate size asymmetry in a binary hard-sphere mixture, as illustrated in Fig. 
EJ It is also the only transition possible for large asymmetries in the case of hard spheres 
provided the system is in the limit of a few large particles dispersed in a sea of small particles 
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(the dynamic arrest of "chancaquilla" limit). 

The second pattern of dynamic arrest is only observed for large enough size disparities 
starting in the regime in which the small particles are a minority. As illustrated in Fig. [21 
here the system may be driven from the region of fully ergodic states (71 = 72 = 00) to 
the region of mixed states (finite 72 and 71 = 00) in which the large particles are arrested 
and the small ones diffuse through the voids left by the large species. In this transition 
72 jumps discontinuously from its infinite value in the fully non-ergodic region to a finite 
value in the region of mixed states, whereas 71 retains its infinite value. Finally, there is 
a transition from this region of mixed states to the region of fully ergodic states in which 
also the small particles are localized (finite values of both, ji and 72). This latter transition 
is characterized by a continuous change in the values of both order parameters 71 and 72, 
although the latter changes continuously from its infinite value inside the region of mixed 
states and right at the transition, to finite values inside the fully non-ergodic region. 

Locating the transitions just described may be done by numerically solving the full self- 
consistent system of equation, as done in Figs. [2] and [3], or by directly determining the non- 
ergodic parameters, which play the role of order parameters of the transition. The latter can 
be done more economically through the solution of the "fixed-point" equations. In section 
IV we derived a strikingly simple and general result for the non-ergodic parameters 71 and 
72, whose use was illustrated with the construction of the phase diagram of the hard-sphere 
binary mixture. We must say that some of the features of the SCGLE results thus determined 



were in fact discovered and discussed in the pioneering work of Refs. [32j, [3J, |35|, |36( based 
on the multi-component version of MCT. Our aim here has been to document what is the 
scenario offered by the SCGLE theory for the same conditions previously studied, and to 
extend the description to new aspects not yet discussed with MCT or with other theoretical 
approaches. In this regard, the main contribution of the present work is the proposal of a 
new theoretical approach that allows the determination of the entire dynamic phase diagram 
of simple model colloidal mixtures. Thus, the second important contribution involves the 
simplest of such models, namely, the hard-sphere binary mixture, for which we outlined the 
entire dynamic arrest phase diagram. In this context we discussed the regime of mild size 
asymmetries, in which we reported the first apparent difference with the results of MCT, 
namely, the fact that the SCGLE theory does not predict the S-shaped transitions predicted 



by Gotze and Voigtmann [37[ on the basis of MCT. For more severe size disparities we also 
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illustrated the main features of the SCGLE predictions, including the precise location of the 
region of mixed states, and its emergence for size disparities beyond a threshold asymmetry 
5 C ~ 0.4. These results do not have yet a MCT counterpart to compare with. 

In this paper we also discussed the fact that other properties besides 71 and 72 could 
be employed as order parameters to classify the various dynamic phases in a mixture. For 
example, one can use the diagonal elements of the non-ergodic parameters associated with 
the self or the collective propagator matrices, which also determine unambiguously the 
various possible states. Although not central to the main subject of this paper, we also 
explained that the use as order parameters of collective correlators, defined as the partial 
intermediate scattering functions divided by their initial value, may lead to a confusing 
description of the hybrid states in which one species is arrested and the other remains 
ergodic. 

In summary, we have presented a theory for dynamic arrest phenomena in colloidal mix- 
tures that offers some advantages over the extension to mixtures of conventional MCT. First, 
it derives from a simpler conceptual formalism, which might allow possible generalizations 
and extensions. Second, the resulting equations are much simpler to solve in practice than 
the corresponding equations of MCT. In fact, the determination of (ergodic-non-ergodic) 
phase diagrams, such as that discussed in the previous section, is considerably simplified 
in the present case through the use of the results for the non-ergodic parameters 7^ in Eq. 
(15.71) . Thus, many important issues regarding the dynamics and the (ergodic-non-ergodic) 
phase behavior of colloidal mixtures may be discussed with the assistance of the SCGLE the- 
ory just presented. The implications of the results presented here in the context of specific 
physical phenomena and experimental conditions will, however, be the subject of separate 
reports. 
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FIG. 1: Time-dependent diffusion coefficients D a (t) as a function of time (upper curves corre- 
sponding to Di(t)), and total intermediate scattering functions Fn(k, t) and F22(k, t) as a function 
of wave- vector of a repulsive Yukawa mixture with z = 0.15, A\ = 10, Ai = 10\/5 for t = to and 
t = 10to- The volume fraction of the more interacting species is kept fixed at 02 = 2.2 x 10~ 4 
and (f>i takes the values 4>i = 0.725 x 10~ 4 (left column), 4>i = 2.2 x 10~ 4 (middle column), and 
4>i = 6.6 x 10 -4 (right column). Solid lines are the SCGLE theoretical results and circles are 
Brownian dynamics data. In the middle column we also include the partial static structure factors 
S aa (k) = F aa (k,t = 0), with the dotted lines being a smooth HNC-fit of the simulation data. 
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FIG. 2: SCGLE results for the collective diffusion propagators ^ a a(k max ,t) at fixed wave- vector 
kmax (the position of the main peak of S aa (k)) for the Yukawa mixture with z = 0.15, §\ = <f>2 = 
2.2 x 10~ 4 calculated using HNC partial static structure factors evaluated at effective coupling 
parameters given by A a = A Q a x, with A® = 10 and A® = 10\/5, for values of x corresponding to the 
reference state in Fig. 1 (x^ = 1.225, curves labeled "1") and to states slightly below (x^ = 4.255, 
labeled "2") and above (x^ = 4.330, labeled "3"), the dynamic arrest transition. The dashed 
plateaus indicate the value of the non-ergodic parameters right at the transition (% = 4.257). 
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FIG. 3: (a) Collective diffusion propagators ^ aa {k max , t) and (b) self diffusion propagators 
^aa{k max , t) at fixed wave- vector k max (the position of the main peak of S aa (k)) for the hard- 
sphere mixture with 5 = 0.3 and x\ = 0.2 for different total packing fraction <\>. The results for 
4> = 0.53 (dotted lines) illustrate the fully ergodic states, whereas those for = 57 (dashed lines) 
illustrate the fully arrested states, in which both species are arrested. The results for (f> = 0.55 
(solid lines) illustrate the hybrid states in which the large particles are arrested while the smaller 
particles continue to diffuse. 
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FIG. 4: Collective diffusion propagators ^ a p{k^t) and diagonal collective diffusion correlators 
&aa(k, t) at a common wave- vector corresponding to the position of the main peak of S22 (k) (see 
the inset, which contains static structure factors at this state) for the hard-sphere mixture with 
5 = 0.3 and X\ = 0.2 for the hybrid state with <\> = 0.55. 
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FIG. 5: Phase diagram of dynamic arrest states of the binary hard-sphere mixture. The three- 
dimensional state space (5, 4>, x\) is spanned by the size disparity parameter 5 = o\j(J2 < 1 5 the total 
volume fraction <j) = 4>i + with cj) a = 7rn Q <7^/6, and the molar fraction of the smaller particles 
x% = n\j{n\ + n,2), where a a and n a are the hard-core diameter and the number concentration 
of particles of species a. This state space is partitioned in three regions, corresponding to fully 
ergodic, mixed, and fully non-ergodic states. In this figure we show the intersection lines of this 
transition surface with the planes of constant 5. The total volume fraction <p has been scaled with 



the volume fraction 



,(m) 



of the ideal glass transition of the mono-disperse hard-sphere system. 



In (a) we illustrate the regime of mild size disparities with the values 5 = 1.0 (horizontal line, 
corresponding to <p = (pg), 0.8, 0.6 (thicker line), and 0.4, where the respective line divides the 
subspace (4>,xi) in the regions of fully ergodic states (below the line) and fully non-ergodic states. 
The four dots linked by the arrow at fixed total volume fraction cfr/cj)^ = 1.0175 correspond to 
the states prepared in the experiment of Williams and van Megen (Ref. 29() for 5 = 0.6. In 



(b) we illustrate the regime of severe size-disparities with the values 5 = 0.3 and 0.2. In this 
regime the region of fully ergodic states also lies below the solid line but a region of mixed states 
appear (shaded region for the case 5 = 0.3). This region is bounded from above by the localization 
transition of the small particles (dashed lines). 
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